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Using 2.92 fb ~ 4 of electron-positron annihilation data collected at y/s = 3.773 GeV with the BESIII detec¬ 
tor, we report the results of a search for the flavor-changing neutral current process D° -A 77 using a double¬ 
tag technique. We find no signal and set an upper limit at 90% confidence level for the branching fraction of 
B(D° -A 77 ) < 3.8 x 10 _G . We also investigate D°-meson decay into two neutral pions, obtaining a branching 
fraction of B(D° -A n°Tv°) = (8.24±0.21(stat.) ±0.30(syst.)) x 10 -4 , the most precise measurement to date 
and consistent with the current world average. 

PACS numbers: 12.60.-i, 13.20.-v, 13.20.Fc, 13.25.Ft 

I. INTRODUCTION the Glashow-Iliopoulos-Maiani mechanism [1], The branch¬ 

ing fraction for D° -A 77 from short-distance contributions. 
In the standard model (SM), the flavor-changing neutral such as an electromagnetic penguin transition, is predicted to 
current (FCNC) decay D° -A 77 is strongly suppressed by be 3 x 10 M- Long-distance contributions due to a vec- 
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tor meson coupling to a photon are expected to enhance the 
branching fraction to the range (1 — 3) x 10 -8 [3, 4]. These 
predictions are orders of magnitude beyond the reach of cur¬ 
rent experiments, but some extensions to the SM can enhance 
FCNC processes by many orders of magnitude. For example, 
in the framework of the minimal supersymmetric SM, gluino 
exchange can increase the branching fraction for the c u 7 
transition to 6 x 10 -6 [5, 6]. 

The previous experimental studies of 79° —>• 77 were per¬ 
formed by the CLEO and BABAR experiments using data sam¬ 
ples collected at the T(4S I ) peak [7, 8]. With an integrated 
luminosity of 470.5 th ' 1 , corresponding to more than 250 
million 79° mesons based on the quoted number of recon¬ 
structed 79° —> TT°n° candidates, its efficiency, and the mea¬ 
sured 75(79° — > 7r°7r°) in Ref. [7], BABAR set an upper limit at 
90% confidence level (CL) on the 79° —> 77 branching frac¬ 
tion of 2.2 x 10 -6 which is the most stringent limit to date. 

In this paper we report a search for 79° — > 77 using 
2.92 ± 0.03 fb _1 of e + e _ annihilation data collected by the 
BESIII detector [9] at ^ = 3.773 GeV in 2010 and 2011. 
There are about 20 million 79° mesons produced [10] from 
•0(3770) decays in this sample. Taking advantage of the 
fact that D-me son production near the •0(3770) resonance is 
solely through 7979, we apply a tagged technique pioneered 
by the MARK III Collaboration [11]. After reconstructing a 
hadronically decaying 79 in an event (the tag), we then search 
for 79-decay candidates of interest in the remainder of the 
event. (Unless otherwise noted, charge conjugate modes are 
implied throughout this paper.) This strategy suppresses back¬ 
ground and provides an absolute normalization for branching 
fraction measurements independent of the integrated luminos¬ 
ity and 7979 production cross section. Therefore, searches for 
D° —> 77 with BESIII at open-charm threshold are uniquely 
clean and provide a valuable complement to studies at the 
T(4S). 

In addition to our primary result, we also report an im¬ 
proved measurement of the branching fraction for the de¬ 
cay D° —» 7r°7r°, which is the dominant background for 
79° —► 77 . Precise measurement of the 79° —> 7 r°7r° branch¬ 
ing fraction can improve understanding of U-spin and SU(3)- 
flavor symmetry breaking effects in 79° decays [12], benefit¬ 
ing theoretical predictions of CP violation in 79 decays [13]. 


II. THE BESIII DETECTOR AND MONTE CARLO 
SIMULATIONS 

The data used in this analysis were collected with the BES- 
III detector operating at the BEPCII Collider. The BESIII de¬ 
tector, which is described in detail elsewhere [14], has a ge¬ 
ometrical acceptance of 93% of 47r and consists of four main 
components. A small-celled, helium-based, multilayer drift 
chamber (MDC) with 43 layers provides momentum resolu¬ 
tion for 1-GeV/c charged particles in a 1-T magnetic field of 
0.5%. Excellent charged particle identification is achieved by 
utilizing the energy loss in the MDC ( dE/dx ). A time-of- 
flight system (TOF) for additional charged particle identifica¬ 
tion is composed of plastic scintillators. The time resolution 


is 80 ps in the barrel and 110 ps in the endcaps, giving 2cr K/ir 
separation for momenta up to about 1 GeV/c. An electromag¬ 
netic calorimeter (EMC) is constructed of 6240 Csl (Tl) crys¬ 
tals arranged in a cylindrical shape (barrel) plus two endcaps. 
For 1.0-GeV photons, the energy resolution is 2.5% in the bar¬ 
rel and 5% in the endcaps. Finally, a muon chamber system 
(MUC) is constructed of resistive plate chambers. These are 
interleaved with the flux-return iron of the superconducting 
magnet. 

Monte Carlo (MC) simulations are used for efficiency 
and background determinations. Events are generated with 
KKMC [15], which incorporates initial-state radiation and the 
spread of the BEPCII beam energy. The generated particles 
are subsequently passed to EVTGEN [16], which simulates 
particle decays based on known branching fractions [17]. To 
realistically mimic our data, we produce a generic MC sample 
including e + e _ —>■ -0(3770) —> DD, continuum hadron 
production (e + e _ —> 7* — > qq, with q = u, d or s), radiative 
returns to the lower cc resonances (e + e _ -A 'yisn (0(3686) 
or J/0)), e + e~ -A t + t~ , and the doubly-radiative Bhabha 
process e + e _ —> e + e~ r y"f. The last component is generated 
with Babayaga [18]. We also generate a signal MC sample 
consisting of e + e _ -A 0(3770) -A 79°79° events in which the 
79° or the 79° decays into a hadronic tag mode or 77, while 
the other 79° or 79° decays without restriction. For all MC 
samples, generated events are processed with GEANT4 [19] 
to simulate the BESIII detector response. 


HI. D° -A 77 ANALYSIS WITH DOUBLE-TAG METHOD 

The 0(3770) resonance is below the threshold for DDn 
production, so the events from e + e _ -A 0(3770) -A DD 
have 79 mesons with energies equal to the beam energy 
(Tt-beam) and known momentum. Thus, to identify 79° can¬ 
didate, we define the two variables AE and Mbc> the beam- 
constrained mass: 

AE = 'y ' Ei — 79beam, 
i 

MBC = ^bea m - I^Z^I 2 ’ 

where E, and p\ are the energies and momenta of the 79° 
decay products in the center-of-mass system of the 0(3770). 
For true 79° candidates, AE will be consistent with zero, and 
Mbc will be consistent with the 79° mass. 

Single tag (ST) candidate events are selected by re¬ 
constructing a 79° in one of the following five hadronic 
final states: 79° -A K + n~, K + n~n 0 , , 

K + 'k~tt + tt~tt°, and 7C + 7r - 7r°7r°, constituting approxi¬ 
mately 37% of all 79° decays [17]. The resolution of Mgg is 
about 2 MeV/c 2 , dominated by the beam-energy spread. The 
AE tag resolutions are about 10 MeV and 15 MeV for final 
states consisting entirely of charged tracks and for those in¬ 
cluding a 7T°, respectively. We search for 79° -a 77 decays in 
these tagged events, thereby highly suppressing backgrounds 
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from QED continuum processes, potential (3770) —► non- 
DD decays, as well as D + D~ decays. The fraction of dou¬ 
ble tag (DT) events, in which the D° is reconstructed as 
D° —> 77, determines the absolute branching fraction for the 
signal mode. 


B{D° -A 77) 


_ 1 y tag,77 _ 

V ]\fi . (A IA \ ' 

1 'tag V c tag,77/ c tag/ 


In this expression i runs over each of the five tag modes, iV Ulg 
and e tag are the ST yield and reconstruction efficiency, and 
iV t ag . 77 and etag ,77 are the yield and efficiency for the DT com¬ 
bination of a hadronic tag and a D° —» 77 decay. 


A. Single-tag selection and yields 


a difference in resolution between data and our MC simula¬ 
tion. Mean and width of the convoluted Gaussian, along with 
the overall normalization, are left free in our nominal fitting 
procedure. The background is described by an ARGUS func¬ 
tion [21], which models combinatorial contributions. In the 
fit, we leave free all parameters of the background function, 
except its endpoint which is fixed at 1.8865 GeV/c 2 . Figure 1 
shows the fits to our tag-candidate samples. Tag yields, given 
in Table I, are obtained by subtracting the fitted background 
estimates from the overall fits in data within the narrow signal 
window Mg® (1.858 < Mg® < 1.874 GeV/c 2 ). The total 
number of tags reconstructed in our data is approximately 2.8 
million. Also shown in Table I are the tagging efficiencies ob¬ 
tained by fitting generic MC Mg^ distributions with the same 
procedure used on data. These ST and DT efficiencies include 
the 7 T° —> 77 branching fraction. 


For each tag mode, D° candidates are reconstructed from 
all possible combinations of final-state particles, according to 
the following selection criteria. Momenta and impact param¬ 
eters of charged tracks are measured by the MDC. Charged 
tracks are required to satisfy |cos 6 *| < 0.93, where 9 is the 
polar angle with respect to the direction of the positron beam, 
and to have a closest approach to the interaction point within 
±10 cm along the beam direction and within 1 cm in the plane 
perpendicular to the beam. Discrimination of charged pions 
from kaons is achieved by combining information about the 
normalized energy deposition ( dE/dx ) in the MDC with the 
flight-time measurement from the TOF. For a positive identi¬ 
fication, the probability of the tt(K) hypothesis is required to 
be larger than that of the I\ (n) hypothesis. 

Electromagnetic showers are reconstructed from clusters of 
energy deposits in the EMC crystals and are required to be in¬ 
consistent with deposition by charged tracks [20]. The energy 
deposited in nearby TOF counters is included to improve the 
reconstruction efficiency and energy resolution. The shower 
energies are required to be greater than 25 MeV for the barrel 
region (|cos0| < 0.80) and greater than 50 MeV for the end- 
caps (0.84 < |cosf?| < 0.92). Showers in the angular range 
between the barrel and endcaps are poorly reconstructed and 
excluded from the analysis. Cluster-timing requirements are 
used to suppress electronic noise and energy deposits unre¬ 
lated to the event. For any tag mode with a n° in the final 
state, photon pairs are used to reconstruct tt° candidates if the 
invariant mass satisfies (115 < ?n 77 < 150) MeV/c 2 . To 
improve resolution and reduce background, we constrain the 
invariant mass of each photon pair to the nominal 7 r° mass. 

For ST modes, we accept D° candidates that satisfy the re¬ 
quirements 1.847 < Mg® < 1.883 GeV/c 2 and \AE^\ < 
0.1 GeV. In events with multiple tag candidates, the one can¬ 
didate per mode with reconstructed energy closest to the beam 
energy is chosen [10]. We extract the ST yield for each tag 
mode and the combined yields of all five modes from fits to 
Mg|, distributions in the samples described above. The sig¬ 
nal shape is derived from the MC simulation which includes 
the effects of beam-energy smearing, initial-state radiation, 
the ±(3770) line shape, and detector resolution. We then 
convolute the line shape with a Gaussian to compensate for 
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FIG. 1 . Fits (solid line) to the M^ c distributions in data (points) for 
the five D° tag modes: (a) K + tt~, (b) K + n~n°, (c) K + iv~ tv + tt ~, 
(d) A' + 7r _ 7r + 7r _ 7r°, and (e) A' + 7r _ 7r°7r°. The gray shaded his¬ 
tograms are arbitrarily scaled generic MC backgrounds. 


B. Double-tag selection and yield 

We select DT candidates by reconstructing D° —> 77 from 
the two most energetic photon candidates that are not used in 
reconstructing the tag mode. The selection criteria for these 
photons are the same as the ones used on the tag side, except 
that we require 0.86 < | cos 6 1 1 < 0.92 for endcap showers 
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TABLE I. Single-tag efficiencies (e[ ag ), tag yields (A tag ) in data, 
double-tag efficiencies (e? a g, 77 ) and their statistical uncertainties. Ef¬ 
ficiencies are determined based on MC simulations. 


modes 

4 (%) <4,77 (%) 

K+ 7T“ 

A' + 7r _ 7r° 

I\ + TV~ 7T + 7T _ 
K + ir~n + Tv~TT° 
K + tv~ 7r°7r° 

66.12 ± 0.04 551800 ± 936 44.8 ± 0.4 

35.06 + 0.02 1097113 ± 1386 24.5 ± 0.1 

39.70 + 0.03 734825 ± 1170 24.7 ± 0.2 

15.32 + 0.04 155899 + 872 9.6 ± 0.1 

15.23 + 0.04 268832 + 976 8.9 ± 0.1 

All Tags 

2808469 ± 2425 


to remove photons landing near the transition region. We re¬ 
quire \AE^\ < 0.10 GeV (1.858 < < 1.874 GeV/c 2 ) 

and \AE" n \ < 0.25 GeV (M^ c > 1.85 GeV/c 2 ) to the 
tag D° candidate and the signal D° candidate, respectively. 
If there are multiple DT candidates, we choose the combi¬ 
nation for which the average of and MgJ (Mbc = 

(Mgg + Mg^,)/2) is closest to the known D° mass [10]. 

For any DT including D° —> K + tt~ , the dominant back¬ 
ground is from the doubly-radiative Bhabha QED process 
e + e _ —> e + e _ 77 , which has a large production cross- 
section. To remove this background, we require the angle be¬ 
tween the direction of the photon candidates and any charged 
tracks to be greater than 10 degrees. This requirement elim¬ 
inates 93% of the QED background. For all tag modes, the 
dominant peaking background in the A E 11 signal region is 
from D° -A 7 r° 7 r°. To remove this background, we implement 
a 7 T° veto. We reject events in which one of the D° -A 77 
final-state photons can be combined with any other photon in 
the event to form a n°. This requirement rejects 82% of the 
D° -a 7 r° 7 r° background and keeps 88 % of the signal events. 
Figure 2 shows the distributions of AE 77 (top) and A E tag 
(bottom) after the above selection criteria are applied, over¬ 
laid with the MC background estimate. 

While we can suppress most of the background with the 
DT method, there remain residual contributions from contin¬ 
uum processes, primarily doubly-radiative Bhabha events for 
Kn tags and e + e _ —>■ qq for other modes. In order to cor¬ 
rectly estimate their sizes, we take a data-driven approach 
by performing an unbinned maximum likelihood fit to the 
two-dimensional distribution of A E 11 versus A /+-. We use 
A E 11 distributions rather than distributions as the back¬ 

ground from non -DD decays is more easily addressed in the 
fit. Also, the background from D° -A n°7r° peaks in Mqq 
at the same place as the signal does, whereas it is shifted 
in AE 1 ' 1 . The fitting ranges are A /s ’ 77 < 0.25 GeV and 
AG lag <0.1 GeV. These wide ranges are chosen to have 
adequate statistics of the continuum backgrounds in our fit. 
The AE 11 resolution is 25 MeV, as determined with signal 
MC. For the signal and the D° -A 7 r° 7 r° background, we ex¬ 
tract probability density functions (PDFs) from MC, where 
the number of D° -A background events is fixed to the 
result of the data-driven method described in Sec. IV. For the 
background from continuum processes, we include a flat com¬ 


ponent in two dimensions, allowing the normalization to float. 
The contribution from D + D~ decays is completely negligi¬ 
ble. We model the background from other D°D° decays with 
a pair of functions. In the A£i tag dimension we use a Crys¬ 
tal Ball Line function (CBL) [22] plus a Gaussian, and in the 
AE 11 dimension, we use a second-order exponential polyno¬ 
mial: 

Y^AE 11 ) = N X g - ( Cl A^ 77 +C2-(AE 77 ) 2 )^ 

In our nominal fitting procedure, we fix the following param¬ 
eters based on MC: the power-law tail parameters of the CBL, 
the coefficients (ci and C 2 ) of the above exponential polyno¬ 
mial, and the mean and the width of the Gaussian function. 
The normalization for the background from all other D°D° 
decays is left free in the fit, as are the mean and width of the 
CBL and the ratio of the areas of the CBL and Gaussian func¬ 
tions. Table I lists the DT signal-reconstruction efficiencies 
for each of the five tag modes. 

As a test to validate the fitting procedure, we fit to 10,000 
sets of pseudo-data (toy MC samples) generated by ran¬ 
domly distributing points based on our generic MC samples 
while taking into account the Poisson distribution with input 
D° -A 77 branching fractions of (0, 5,10) x 10 -6 . The av¬ 
erage branching fractions measured with these samples are 
(0.3 ± 1.2,5.0 ± 2.4,10.0 ± 3.1) x 10 -6 , respectively, where 
the quoted uncertainties are the root-mean-squares of the dis¬ 
tributions. 

Figure 2 shows projections of the fit to the DT data sample 
onto AE 11 (top) and A E tag (bottom). We also overlay back¬ 
ground distributions predicted by the MC simulations. The 
fit yields Mag ,77 = (— 1.04 3)’ demonstrating that there is 
no signal for D° -A 77 in our data. This corresponds to 
B(D° -A 77 ) = (— 0 . 6 ^ 3 ) x 10 ~ 6 where the uncertainties 
are statistical only. 


IV. SIZE OF D° -a 7r°7r° BACKGROUND 


To estimate the contribution of background from D° -a 
7 r n 7 T° events to our selection, we make a second DT measure¬ 
ment with the same sample used in searching for D° -A 77 . 
Within these tagged events, we reconstruct D° -A 7 r° 7 r° with 
the 7 r° candidates that are not used in reconstructing the tag 
modes. The selection criteria for these 7 r° candidates are 
the same as those used in reconstructing the tags. We se¬ 
lect the pair of 7 r°s that gives the smallest \AE n 77 | and 


extract the DT yield by fitting to MgJ 7 , while requiring 
—0.070 < A E n °*° < +0.075 GeV. In this fit, a double- 
Gaussian function is used to represent the shape for 

the D° -A 7T 0 7r° decays, while the D°D° MC shape describes 
the background. 


Figure 3 shows the fit to the Mg C 
M,T>° < 1.886 GeV/c 2 , which yields N ob 0 s 0 = 1036 ± 35 


distribution in 1.840 < 

T°7r° 

J BC 

events for D° -a Tr°n°. Thus the yield in our data sample 
of D° -A 7r°7r° with a D° decaying into one of the five tag 
modes is = IV^o/e^ 770 , where = 6.08% is 

the DT efficiency for D° -a 7r°7r° as determined with MC. 
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FIG. 2. Fit to the DT sample in data (points), projected onto AE 77 
(a) and AE tag (b). The dashed lines show the overall fits, while 
the dotted histograms represent the estimated background contri¬ 
bution from D° —> 7 r° 7 r°. The solid line superimposed on the 
AE 77 projection indicates the expected signal for B(D° — > 77 ) = 
10 x 10 -6 . Also overlaid are the overall MC-estimated backgrounds 
(gray shaded histograms) and the background component from non- 
HD processes (diagonally hatched histograms). 


The expected 7 r° 7 r° contribution to our 77 candidates can be 
then obtained as 

,77 

,, rexpected _ ^produced 77 _ xrobs t 7 r° 7 r° 

7T°7T° 7T 0 7T 0 7r°7T 0 7T° 7T° ^0^0 

e DT 

where = 0.11% is the efficiency for D° -7 7 r° 7 r° to 

be counted as D° —> 77 . The efficiencies e^o^o and 
include the reconstruction efficiencies for the tag sides as well 
as the branching fractions, although these cancel in the ratio. 

We consider the following sources of systematic uncer¬ 
tainty in determining the D° —> 7 r° 7 r° contamination: ir° re¬ 
construction (1.5%), photon reconstruction (2.0%), binning of 
^b°g° (0-1%). fit range ( 0 . 1 %), background shape (0.5%), 
signal shape (1.7%), and the A E* n requirement (0.6%). 


FIG. 3. Fit to the distribution in data (points) for D° —> 

7T 0 7r° DT candidates. The solid line is the total fitted result, while 
the dotted and dashed lines are the background and signal compo¬ 
nents of the fit, respectively. The diagonally shaded histogram is the 
background determined with MC. 


Combining statistical and systematic uncertainties, we esti¬ 
mate the number of D° —» events among the D° —► 77 

candidates to be 18 events with a relative uncertainty of 4.6%, 
spread across the A E 11 fit range. 


V. SYSTEMATIC UNCERTAINTIES FOR D° -a 77 
ANALYSIS 


MC studies demonstrate that D-decay measurements based 
on DT-to-ST ratios benefit from cancellation of most of the 
systematic uncertainties of tag reconstruction. The overall 
systematic uncertainty in our measurement is therefore domi¬ 
nated by other effects. The systematic uncertainties that are 
independent of our signal-fitting procedure are that associ¬ 
ated with detection of the two photons, which is estimated 
by studying the reconstruction efficiency of a daughter pho¬ 
ton from 7 T° decay in a DT D° -A Kg 7 r° sample (2.0%); 
the signal-side requirement, which is estimated from 

the AE 7r w distribution of the DT D° -A 7 r° 7 r° sample and 
by observing the stability of the B{D° -A 7 r° 7 r°) while vary¬ 
ing the selected range of Mg C 7r (3.1%). The systematic un¬ 
certainties in ST yields (1.0%) are estimated first for individ¬ 
ual tag modes, and then combined in quadrature with weights 
based on the observed tag yields (iVA). The sources for the 
uncertainties of ST yields we consider are the choice of fit 
range, assumed signal parametrization, and the signal 

window. Combined in quadrature, these total 3.8%. 

We also consider six possible sources of systematic 
effects due to our fitting procedure. (i) Fits are re¬ 
done with all possible combinations of fitting ranges: 
— (0.12,0.10,0.08)<AE tag <+(0.08,0.10,0.12) GeV and 
-(0.30,0.25,0.20)<AET'T<+(0.20,0.25,0.30) GeV. (ii) 
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The MC-based analytic form of the D°D° background 
shape (excluding the D° —► 7r°7r n contribution) is var¬ 
ied by changing the input branching fractions for D° —► 
tt° ij/r]ri/K^r//K^n° by ± 1 < 7 pdg [17]. (iii) The flat non -DD 
background shape is replaced with a shape that is linear in the 
A E 11 dimension, (iv) The fixed size of the background from 
D° —► 7r°7r° is varied by ±4.6%. (v) The fixed shape of the 
background from D° —> n°n° is studied by comparing A E 
distributions of DT events from D° —> / K^n 0 / Ktttt 0 

between data and MC simulations in which we intention¬ 
ally ignore the lower-energy photon from each 7r° decay to 
mimic our background. We conclude that we do not need to 
assign additional systematic uncertainty due to the assumed 
D° —> 7r°7r° background shape in the fit, except to give an 
extra Gaussian smearing of a = 5 MeV in the AP ulg dimen¬ 
sion. (vi) The fixed signal shape is studied based on the DT 
D° —► TT°n° sample in which we study distributions of its 
A P Ulg and A E n ~ for four cases by requiring that one of the 
two photons from each of the two ir° to have at least 0.5, 0.6, 
0.7, and 0.8 GeV to mimic our signal photon energies. From 
all four cases, we find that we need an extra Gaussian smear¬ 
ing of cr = 16 MeV and a shift by a factor of 1.0025 in the 
A E 11 dimension as well as an extra smearing of a 5 MeV 
in the AP Ulg dimension. 

Table II summarizes systematic uncertainties that are inde¬ 
pendent of our fitting procedure, as well as systematic vari¬ 
ations that we consider to estimate uncertainties due to the 
fitting procedure. In the next section, we describe how we 
combine these systematic uncertainties into our measurement. 

TABLE II. Systematic uncertainties and variations for D° —> 77 
analysis. 


Uncertainties independent of fitting procedure 


Source 

Relative uncertainty (%) 

Photon reconstruction 

2.0 

requirement 

3.1 

ST D° yields 

1.0 

Total 

3.8 


Systematic variations due to fitting procedure 


Source 

Variations 

Fit range (GeV) 

±0.02 in P tag and ±0.05 in P 77 

D° —> 7 r° 7 r° norm. 

±4.6% 

D° 7 r° 7 r° shape 

Smear in AP tag 

D°D° bkg shape 

ASinput[P° -A (rin 0 /r/ri/K^n 0 /K^ri)) 

Non-D°Z ) 0 bkg shape 

Flat vs Linear 

Signal shape 

Smear in AP tag and AP 77 , shift in P 77 


VI. THE RESULT FOR D° -> 77 

Since we do not observe a signal, we set an upper limit on 
the branching fraction for D° —» 77 . We first obtain a smooth 


background-only PDF shape from the sample via the kernel 
estimation method [23], This is done by utilizing the RooFit 
class [24] RooNDKeysPdf [25], We then generate 2.2 million 
toy MC samples by randomly distributing points according to 
the PDF shape, while taking into account the Poisson distri¬ 
bution. We fit to each of these toy samples while randomly 
making systematic variations in the fitting procedure, as de¬ 
scribed in the previous section. We also simultaneously smear 
each of the fitted branching fractions with a Gaussian whose 
width (3.8%) corresponds to the total systematic uncertainty 
that is not associated with the fitting procedure. 

Figure 4 shows an accumulation of the resulting branching 
fractions for D° —> 77 . The shaded region represents 90% 
of its physical region, which we use to set our 90% CL upper 
limit of B(D° -A 77 ) < 3.8 x 10 -6 . If the systematic uncer¬ 
tainty were ignored in setting this limit it would be reduced 
by 0.1 x 10 6 . The expected measurement of branching frac¬ 
tion from these toy experiments is x 10 -6 , where 

the quoted uncertainties correspond to 68 % of the areas un¬ 
der the curves in Fig. 4. The mean value of the accumulated 
branching fractions is consistent with the value of the branch¬ 
ing fraction from the nominal fit to data at 0 . 6 <r level. 
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FIG. 4. Accumulated branching fraction distribution based on toy 
MC samples generated from the data-driven PDF. (See the text for 
details.) The shaded region represents 90% of the physical region. 



VII. IMPROVED MEASUREMENT OF D° ->• 7r°7r° 
BRANCHING FRACTION 


As a byproduct of this analysis we also measure the branch¬ 
ing fraction of D° —> 7 r° 7 r° using the same data sample. Since 
the produced D°D° pairs in our sample necessarily have op¬ 
posite CP eigenvalues [ 20 ], the effective branching fraction 
for the CP-even final state 7 r° 7 r° is altered when it is mea¬ 
sured in events tagged with a CP-mixed state such as D° —> 
K+tt~ [26]. To avoid this complication and to improve the 
statistics, instead of a DT technique, we reconstruct only one 
D° or D° decay in the ip(3770) -A D°D° process. The ob- 
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served yield is normalized to the total number of the D°D° 
pairs, which can be obtained as N d oqo = £ x a(e + e~ —► 
■0(3770) — > D°D°), using the integrated luminosity C of 
our sample [9] and the previously measured cross section 
a{e+e~ -> D°D°) = (3.607±0.017(stat.)±0.056(syst.)) nb 
[10]. The branching fraction for D° —> 7 r° 7 r° can be calcu¬ 
lated as 


B(D° -+ 7r°7T 0 ) 


^7r°7r° 

e-n-OTT 0 • 2-ZVpopo 


where N^o^o is the observed number of D° -+ tt°tt 0 decays 
and e^o^o is the selection efficiency determined with MC. 

The reconstruction of 7 r° candidates is the same as those 
in the ST modes described in Sec. IIIA. We choose a pair of 
reconstructed 7 r°s that give the smallest | A E v 7r |, and require 
—0.06 < A E* 0 * 0 < +0.03 GeV. The resolution of AE n ° n ° 
is about 20 MeV. Then we extract the signal yield from a fit 
to 5 C 7r . The efficiency is determined to be e n o v o = 36% 
from MC simulations. 

Figure 5 shows a fit to the distribution in 1.8400 < 

Mg,+ < 1-8865 GeV/c 2 . We use a double-Gaussian func¬ 
tion to describe the signal shape, which is shown as a dot¬ 
ted line, and the background shape is described by an AR¬ 
GUS background function [21]. From this fit, which yields 
X 2 /d.o.f. = 91.8/85, we obtain N n o n o = 6277 ± 156 events. 
In Fig. 5, we also overlay the backgrounds that are estimated 
by the MC simulations (gray shaded histogram). 



FIG. 5. Fit to the distribution in data for D° n°n° candi¬ 

dates (points). The shaded histogram is the background predicted by 
MC. The solid and dashed curves are the total fit and the background 
component, respectively, and the dotted curve shows the signal. 


uncertainty due to n° reconstruction is estimated with a DT 
D° -+ A' _ 7 r + 7 r° sample, (ii) Histogram binning scheme is 
varied, (iii) Narrower (1.8450 < M^ ( f < 1.8820 GeV/c 2 ) 
and broader (1.8350 < Mgg r ° < 1.8865 GeV/c 2 ) fit ranges 
are tried, (iv) Narrower (-0.055 < A E*°”° < 0.025 GeV) 
and broader (—0.065 < AH 7r ° 7r ° < 0.035 GeV) requirements 
are applied, (v) Instead of using the ARGUS function [21], a 
MC-based background shape is used, (vi) To assess a posible 
bias due to the signal line shape, we fix the all shape parame¬ 
ters of the double Gaussians based on the shape extracted from 
the DT D° —» K~tt + tt° sample, (vii) The uncertainty of the 
determination of N D oqo is determined based on Refs. [9, 10]. 
The resultant relative uncertainties are shown in Table III. 


TABLE III. Systematic uncertainties for D° —> 7r°7r 0 analysis. 


Source 

Relative uncertainty (%) 

7r° reconstruction 

1.5 

Histogram binning 

0.1 

Fit range 

2.4 

A E™ n requirement 

0.6 

Background shape 

0.2 

Signal shape 

0.9 

Nd°d° 

1.9 

Total 

3.6 


VIII. CONCLUSIONS 

Using 2.92 flu 1 of e + e _ annihilation data collected at 
= 3.773 GeV with the BESIII detector, we have searched 
for the FCNC decay D° -+ 77 and observe no significant sig¬ 
nal. We set an upper limit B(D° —> 77 ) < 3.8 x 10 -6 at 
the 90% CL, which is consistent with the upper limit previ¬ 
ously set by the BABAR Collaboration [7] and with the SM 
prediction. Ours is the first experimental study of this decay 
using data at open-charm threshold. Employing the DT tech¬ 
nique, we are able to suppress the backgrounds from non-Hi) 
decays effectively. Our analysis also shows that the peaking 
background from D° -+ 7 r° 7 r° can be reliably estimated with 
a data-driven method. 

We have also measured the branching fraction for D° —> 
7 r° 7 r° to be (8.24 ± 0.21(stat.) ± 0.30(syst.)) x 10 -4 which is 
consistent with the previous measurements [27] and the most 
precise to date. 


From the fitted signal yields (A+n_o) and reconstruction ef¬ 
ficiency (e^o^o), we obtain 

B{D° -+ ttV) = (8.24 ± 0 . 21 (stat.) ± 0.30(syst.)) x 10" 4 . 

The quoted total systematic uncertainty (3.6%) is the quadra¬ 
ture sum of the following seven sources of uncertainty, (i) The 
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